POPULAR SUMMARY 


Mean ages of stratospheric air derived from in situ observations of C0 2 , CH 4 , and N 2 0 

Carbon dioxide concentrations at the Earth's surface show a long-term increase due to 
fossil fuel combustion. In addition, carbon dioxide varies seasonally because it is removed 
from the atmosphere during the growing season by photosynthetic plants. Seasona 
differences depend on the type and amount of vegetation and on other factors, such as 
temperature, rainfall, and the length of the day. The seasonal cycle and long-term increase 
propagate from the surface to high altitudes, eventually reaching the stratosphere. 
Measurements of carbon dioxide can therefore be used to study the rate at which air moves 

throughout the atmosphere. 

Data from instruments flown on board the NASA ER-2 aircraft and high-altitude 
balloons have been used to create a map of the average age of stratospheric air, whic 
provides information about how long pollutants will remain in the atmosphere The data 
were collected from 1992 through 1998 with nearly pole-to-pole coverage, including all 
seasons except Southern Hemisphere summer. The resulting comprehensive recor o 
stratospheric carbon dioxide concentrations is unprecedented and extremely precise. e 
instruments can measure a change in the carbon dioxide concentration of one part in ten 
million (or 0.00001 percent). The distribution of the average age of stratospheric air wi 
serve as a critical test for models used to understand the response of the stratospheric ozone 
layer to natural and anthropogenic perturbations, such as volcanic eruptions, the buildup an 
subsequent gradual decrease of chlorofluorocarbons (CFCs), or emissions from high-flying 
aircraft like the Concorde. This dataset provides a benchmark for measuring future changes 
in the stratospheric circulation that might result from climate change or other processes. 


A. E. Andrews, K. A. Boering, B. C. Daube, S. C. Wofsy, M. Loewenstein, H. Jost, 

J. R. Podolske, C. R. Webster, R. L. Herman, D. C. Scott, G. J. Flesch, E. J. Moyer, 

J. W. Elkins, G. S. Dutton, D. F. Hurst, F. L. Moore, E. A. Ray, P. A. Romashkin, and S. E 

Strahan 



Mean ages of stratospheric air derived from in situ 
observations of CO2, CH4, and N2O 

A. E. Andrews , 1 ’ 2 K. A. Boering , 2 B. C. Daube , 1 S. C. Wofsy, M. Loewenstein, 
H. Jost , 4 ’ 5 J. R. Podolske , 4 C. R. Webster , 6 R. L. Herman , 6 D. C. Scott , 6 
G. J. Flesch , 6 E. J. Moyer , 6 ' 7 J. W. Elkins , 8 G. S. Dutton , 8 D. Hurst , 8 
F. L. Moore , 8 E. A. Ray , 8 P. A. Romashkin , 8 and S. E. Strahan 9 


Abstract. Accurate mean ages for stratospheric air have been derived from a spatially and 
temporally comprehensive set of in situ observations of C0 2 , CH 4 , and N 2 0 obtained from 
1992 to 1998 from the NASA ER-2 aircraft and balloon flights. Errors associated with the 
tropospheric C0 2 seasonal cycle and interannual variations in the C0 2 growth rate are < 0.5 
year throughout the stratosphere and < 0.3 year for air older than 2 years (N 2 0 < 275 ppbv), 
indicating that the age spectra are broad enough to attenuate these influences over the time pe- 
riod covered by these observations. The distribution of mean age with latitude and altitude 
provides detailed, quantitative information about the general circulation of the stratosphere. At 
20 km, sharp meridional gradients in the mean age are observed across the subtropics. Between 
20 and 30 km, the average difference in mean age between the tropics and midlatitudes is 
~2 years, with slightly smaller differences at higher and lower altitudes. The mean age in the 
midlatitude middle stratosphere (-25-32 km) is relatively constant with respect to altitude at 
5 ± 0.5 years. Comparison with earlier balloon observations of C0 2 dating back to the 1970s 
indicates that the mean age of air in this region has remained within ±1 year of its current value 
over the last 25 years. A climatology of mean age is derived from the observed compact rela- 
tionship between mean age and N 2 0. These characteristics of the distribution of mean age in 
the stratosphere will serve as critically needed diagnostics for models of stratospheric trans- 
port. 


1. Introduction 

Transport of trace species into and within the stratosphere 
must be accurately represented in models of the stratospheric 
ozone layer in order to reliably simulate the atmosphere’s re- 
sponse to natural and anthropogenic perturbations. Air enters 
the stratosphere primarily through the tropical tropopause, rises 
to higher altitudes, then moves poleward and descends at mid- 
dle and high latitudes, eventually returning to the troposphere 
[e.g., Brewer , 1949; Dobson, 1956; Holton et ai , 1995]. 
However, many details of the circulation, such as the altitude 
dependence of horizontal transport rates and the relative impor- 
tance of advective and diffusive processes, are still not well 
quantified, and inadequate representation of constituent trans- 
port is perhaps the largest single source of uncertainty in cur- 
rent models of the ozone layer [Kawa etal, 1999]. 

The “mean age” of stratospheric air has become an impor- 
tant diagnostic for evaluating transport simulations in models 
of the stratosphere. A stratospheric air parcel may be consid- 
ered to be comprised of infinitesimal fluid elements with di- 
verse transport histories since crossing the tropical tropopause 
[ Kida , 1983; Hall and Plumb , 1994]. The statistical distribu- 
tion of transit times for the fluid elements in the parcel is called 
the “age spectrum,” equivalent to the Green’s function between 



the boundary at the tropical tropopause and the interior of the 
stratosphere, and the mean age corresponds to the first moment 
of this distribution [Hall and Plumb , 1994]. Thus the mean 
age of the parcel is the average transit time from the tropical 
tropopause, calculated over the ensemble of fluid elements 
comprising the parcel. 

For an inert tracer with a linearly increasing source, the 
mean age at a given point in the stratosphere would be equiva- 
lent to the elapsed time between the occurrence of a particular 
mixing ratio at the tropical tropopause and the occurrence of 
the same mixing ratio at the point in the stratosphere [Hall and 
Plumb , 1994]. Several studies have shown that mean ages can 
be estimated from observations of nearly linearly increasing 
tracers such as CO 2 [e.g., Bischof et al, 1980, 1985; Schmidt 
and Khedim , 1991; Nakazawa et al , 1995; Boering et al, 
1996; Strunk et al. , 2000] or SF 6 [e.g., Maiss et al. , 1996; 
Harnisch et al. , 1996; Volk et al, 1997; Strunk et al, 2000]. 
The recent NASA Models and Measurements II (MMII) study 
revealed a large discrepancy between mean ages predicted by 
the majority of two-dimensional (2-D) and three-dimensional 
(3-D) models of the atmosphere and those derived from obser- 
vations of SFg and CO 2 , with most models significantly under- 
estimating mean ages compared to observations [Park et al, 
1999; Hall et al , 1999]. Several studies have investigated 
factors controlling mean age in these models and the implica- 
tions for predicting the stratospheric response to natural or 
anthropogenic perturbations [e.g., Bacmeister et al , 1998; 
Fleming et al, 1999; Li and Waugh , 1999; Hall and Waugh, 
2000; Fleming et al, 2001; Jones et al, 2001]. However, 
uncertainties regarding the nature and significance of the dis- 
crepancy between modeled and tracer-based mean ages re- 
main. 

Boering et al [1996] noted that mean age is mathematically 
equivalent to a conserved tracer with a uniform stratospheric 
source of 1 yr/yr and may therefore be viewed as an approxi- 
mate surrogate for exhaust that would be emitted by strato- 
spheric aircraft. Although emissions from any future strato- 
spheric aircraft are likely to be concentrated in the lower strato- 
sphere along certain flight corridors, rather than to be emitted 
homogeneously throughout the stratosphere, mean age is more 
representative of aircraft exhaust than most other common 
tracers of stratospheric transport (e.g., N 2 0, CH 4 , and the 
CFCs), which are characterized by tropospheric sources and 
stratospheric sinks. Hall and Waugh [2000] used a simple 
model to further explore the relationship between mean age 
and the residence time of stratospheric pollutants and found 
that mean age and residence time were correlated, although the 
correlation was imperfect. In particular, the stratospheric resi- 
dence time was quite sensitive to varying the height of the mid- 
latitude tropopause, while mean age was largely unaffected. 

The results of MMII showed that models with systemati- 
cally young mean ages tended to have more N 2 0 and less NO v 
and Cl v than those with older mean ages, raising concern that 
models which underestimate mean age may also underestimate 
the impact of chemical perturbations to the stratosphere, such 



as an increase in NO v resulting from a fleet of stratospheric 
aircraft [Park et ai , 1999; Hall et ai, 1999). Li and Waugh 
[1999] obtained similar results using a 2-D model to explore 
the sensitivity of mean age and long-lived tracers to variations 
in the transport parameters. 

Recently, Fleming et al. [2001] compared results from the 
1995 and 1999 versions of the NASA Goddard Space Flight 
Center 2-D model and showed that mean ages in the 1999 
version agree much better with mean ages derived from obser- 
vations. The total ozone simulations for both versions were in 
reasonably good agreement with each other and with data from 
the Total Ozone Mapping Spectrometer, but the modeled ozone 
loss resulting from a perturbation due to exhaust from strato- 
spheric aircraft was a factor of 2.4 larger in the 1999 version, 
illustrating that simulations of inert tracers can reveal impor- 
tant model transport deficiencies that are not apparent in simu- 
lations of chemically active species. 

In this paper, we present stratospheric mean ages derived 
from the extensive set of in situ observations of stratospheric 
C0 2 , CH 4 , and N 2 0 obtained from 1992 to 1998 using the 
NASA ER-2 aircraft and high-altitude balloons. This work is 
intended to provide accurate nearly global distributions of 
stratospheric mean age. Errors in the derived mean ages asso- 
ciated with the seasonal cycle in C0 2 mixing ratios in air enter- 
ing the stratosphere and with interannual variations in the tro- 
pospheric C0 2 growth rate are analyzed. The associated errors 
are shown to be insignificant for air with N : 0 mixing ratios 
< 275 ppbv, corresponding to mean ages > -1.9 years. For air 
with N 2 0 > 275 ppbv, for which errors in mean age due to 
seasonal and interannual variations in C0 2 are significant, we 
use the empirical age spectra derived from C0 2 observations 
by Andrews et ai [1999, 2001] to obtain accurate values for 
mean age. These mean ages from C0 2 agree with mean ages 
derived from simultaneous measurements of SF 6 to within 20% 
and are also consistent with mean ages derived from earlier 
measurements of both C0 2 and SF 6 in whole air samples ob- 
tained using balloon-borne cryosamplers. Finally, a compari- 
son of published C0 2 mixing ratios measured in the midlati- 
tude middle stratosphere since the mid-1970s shows that the 
mean age in this region has been relatively constant over the 
past 25 years. 

2. Observations 

The measurements used in this analysis were obtained on 
flights of the NASA ER-2 aircraft and high-altitude balloons. 
The ER-2 measurement campaigns occurred from 1992 to 
1998 and included the Stratospheric Photochemistry and Dy- 
namics Expedition (SPADE), the Airborne Southern Hemi- 
sphere Ozone Experiment/Measurements for Assessing the 
Effects of Stratospheric Aircraft (ASHOE/MAESA), Strato- 
spheric Tracers of Atmospheric Transport (STRAT), Photo- 
chemistry of Ozone Loss in the Arctic Region in Summer 
(POLARIS), and a short test flight series in November 1998. 
The balloon flights were conducted during the Observations of 


the Middle Stratosphere (OMS) experiment, which consisted 
of deployments from Fort Sumner, New Mexico (35°N); 
Juazeiro do Norte, Ceara, Brazil (7°S); and Fairbanks, Alaska 
(65 °N). The maximum altitude of the ER-2 aircraft is ~21 km, 
while the balloons typically reached altitudes of -32 km. De- 
tailed descriptions of the dates and latitude ranges for each 
deployment are given in Table 1 for the ER-2 deployments and 
in Table 2 for the balloon flights. This set of data presents a 
uniquely comprehensive picture of stratospheric composition in 
time and space, allowing the analysis of mean age to account 
for spatial, seasonal, and interannual variations that could not 
be addressed using previously available data. 

Mixing ratios of C0 2 were measured by nondispersive infra- 
red absorption by two separate instruments, one on each plat- 
form. The design of the balloon C0 2 analyzer was based on 
that of the ER-2 instrument, with modifications allowing it to 
operate at much higher altitudes (B. C. Daube et al., A high- 
precision fast-response airborne C0 2 analyzer for in situ sam- 
pling from the surface to the middle stratosphere, submitted to 
Journal of Atmospheric Oceanic Technology , 2001) (hereinaf- 
ter referred to as (Daube et ah, submitted manuscript, 2001)). 
Note that two of the balloon flights, on September 21, 1996 
and June 30, 1997, were nearly spatially and temporally coin- 
cident with ER-2 flights. Data from these flights indicated 
excellent agreement between the balloon and ER-2 instru- 
ments. The June 1997 flight showed that the balloon and ER-2 
C0 2 instruments agreed to within 0.1 ppmv. 

Both C0 2 instruments are calibrated in flight, giving a 
long-term precision better than 0.1 ppmv from 1992 to 1998 
for all deployments except the September and December 1996 
deployments of the STRAT experiment, when degradation of a 
component in the detector led to slightly reduced precision of 
-0.15 ppmv. The calibration standards are traceable to stan- 
dards held by the Scripps Institute of Oceanography (SIO) and 
to a second set of standards calibrated by the National Oceanic 
and Atmospheric Administration/Climate Monitoring and 
Diagnostic Laboratory (NOAA/CMDL). Stratospheric data 
have been reported on the SIO X95 mole fraction scale. We 
have participated in two round robin intercomparisons organ- 
ized by the World Meteorological Organization, so the data 
can be related to that from the other participating institutions 
[Francey, 1997]. The NOAA/CMDL scale differs by up to 
0.2 ppmv from the SIO X95 scale for the relevant range of C0 2 
mixing ratios, with the SIO X95 scale having higher values. 
However, the maximum difference between the SIO X95 and 
NOAA/CMDL scales is equivalent to a 0.1 year difference in 
mean age and will be neglected for this analysis. 

In the stratosphere, C 0 2 is conserved, apart from a small 
source from CH4 oxidation. Measurements of CH4 are avail- 
able on both platforms and are used to correct for stratospheric 
production of C 0 2 , by assuming that the difference between 
observed CH4 and the tropospheric value (-1.725 ppmv for this 
time period) corresponds to the amount of C 0 2 produced. CH4 
was measured by the Aircraft Laser Infrared Absorption Spec- 
trometer (ALIAS) on the ER-2 with a stated accuracy of 5% 



[Webster et ai , 1994]. Since CH 4 data are not available for all 
flights, we use an average relationship between ALIAS CH4 
and N 2 0 measurements from the Airborne Tunable Laser Ab- 
sorption Spectrometer (ATLAS) on the ER-2 (see below): 

CH 4 = 0.56 + 0.002287 N 2 0 + 4.6 x 1 0 6 N 2 0 2 , ( 1 ) 

where CR* is given in ppmv and N 2 0 is in ppbv. For flights 
where CR, data are available, differences between measured 
CH 4 and the value predicted by ( 1 ) are negligible in terms of 
the C0 2 source correction. Therefore for simplicity we use (1) 
and measured N 2 0 to correct for the stratospheric production of 
C0 2 from CH 4 for all ER-2 flights. Equation (1) is robust over 
the 6 year time period of the measurements used in this analy- 
sis despite increases in the tropospheric mixing ratios of both 
N 2 0 and CH4 during this time, because the growth rates of 
both species from 1992 to 1998 are relatively small. Note that 
for other aspects of this analysis, N 2 0 mixing ratios have been 
corrected for growth, as described below. 

On the balloon flights, CH4 was measured by the Airborne 
Laser Infrared Absorption Spectrometer (ALIAS-II) [Scott et 
al. , 1999] with a measurement accuracy of 5%. CH* data from 
ALIAS-II are not available for the November 1 1, 1997, flight 
from Juazeiro do Norte, Brazil (7°S). For this flight, methane 
mixing ratios were estimated using N 2 0 data from the Argus 
instrument (see below) and the tropical CH 4 iN 2 0 relationship 
determined from the Atmospheric Trace Molecule Spectros- 
copy Experiment (ATMOS) [Michelsen et ai, 1998]. Several 
filaments of air which appeared to have been recently trans- 
ported from midlatitudes were observed during this flight [Jost 
et ai , 1998], for which the midlatitude CH 4 :N 2 0 relationship 
from ATMOS may provide more accurate estimates of CH4. 
However, the difference for the C0 2 correction between the 
tropical and midlatitude CH^N-O correlations from ATMOS 
is negligible for this analysis (<0.15 ppmv) for the relevant 
range of N 2 0 mixing ratios. 

N 2 0 was measured on the ER-2 by the ATLAS instrument 
as described by Podolske and Loewenstein [1993] and on the 
OMS balloon gondola by ALIAS-II and Argus [Jost et ai , 
1998]. The measurement accuracy for ATLAS and Argus is 
2.5%, while that for ALIAS and ALIAS-II is 5%. We cor- 
rected for growth in the tropospheric mixing ratio of N z O from 
1992 to 1997 by scaling data from each flight by a factor of 
1 . 0020 (l997t) , where t is the flight date in years, following 
Strahan et ai [1999]. This factor corresponds to an average 
trend of 0 . 2 % per year. 

SF 6 and CFC-1 1 were measured by gas chromatography by 
the four channel Airborne Chromatograph for Atmospheric 
Trace Species (ACATS-IV) on the ER-2 [Elkins et ai , 1996; 
Romashkin et ai , 2001] and by the Lightweight Airborne 
Chromatography Experiment during the balloon flights [Ray et 
ai, 1999]. The accuracy of the measurements typically ranges 
from 1 to 4% depending on the species and the flight [Ray et 
ai, 1999]. 
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3. Analysis 

Detailed knowledge of the stratospheric boundary condition 
for C0 2 (i.e., the C0 2 mixing ratio in air entering the strato- 
sphere) is required to calculate mean age from observations of 
stratospheric C0 2 mixing ratios. Newly stratospheric air en- 
countered by the ER-2 was identified by near-tropospheric N 2 0 
mixing ratios (-310 ppbv), high CO (30 to 50 ppbv), low H 2 0 
« -6 ppmv), and potential temperature > 390 K. (In this 
analysis, we consider only the “stratospheric overworld,” as 
discussed by Holton et ai [1995], where isentropes are con- 
tained completely within the stratosphere.) Air with these 
characteristics was frequently sampled by the ER-2 near the 
tropical tropopause and at higher latitudes in filaments of air 
recently transported from the tropics. A continuous representa- 
tion for the observed boundary condition was derived (Figure 
la), assuming that the seasonal cycle is invariant and that the 
long-term trend can be represented by the deseasonalized 
(12-month running mean) average of surface data from 
Mauna Loa (19°N) and Samoa (14°S) delayed by 2 months, as 
shown by Boering et al. [1996]; see Andrews et al. [1999] for 
additional details. Note that the deseasonalized average of 
data from Mauna Loa and Samoa represents the global surface 
mean given by Conway et ai [1994] to a very good 
approximation. 

The elapsed time since the C0 2 mixing ratio observed in an 
air parcel (corrected for CH4 oxidation) last occurred at the 
tropical tropopause provides a simple estimate of the mean age 
F For clarity, estimates of the mean age derived in this man- 
ner will be called C0 2 lag times r LAG . Initially, we will ignore 
seasonal and interannual variations in C0 2 and calculate /Lag 
using a linear fit to the long-term trend in C0 2 mixing ratios 
from 1980 to 1998: 

r LA G = t - (1749.88 + 0.68241 (C0 2 - (1.725 - CH*))). (2) 

Here, t is the date of the observation in years, C0 2 and CH4 are 
in units of ppmv, and the mean time to add 1 ppmv C0 2 is 
0.68241 year (corresponding to a growth rate of 
1.47 ppmv yr' 1 ). Note that since we do not have direct 
observations of the stratospheric boundary condition prior to 
November 1992, we must assume that the relationship 
between the stratospheric boundary condition and surface 
station data has not changed significantly since 1980. 
Published observations of stratospheric C0 2 from balloons 
since the 1970s indicate that this assumption is valid (see 
sedfditad.fefi^.of increase of C0 2 were constant at 1.47 ppmv 
yr' 1 and there were no seasonal cycle, r ^ a g would exactly equal 
the mean age. However, the rate of increase of C0 2 has varied 
considerably over the last two decades [Keeling et ai , 1995], 
as can be seen by comparing the long-term trend derived from 
surface data with the linear fit used to calculate r LAG (Fig- 
ure la). In fact, the growth rate of the deseasonalized 
(12-month running average) continuous stratospheric boundary 
condition varies from 0.5 to >3 ppmv yr 1 over this time period 
(Figure lb). Furthermore, there is a strong seasonal cycle in 
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C0 2 mixing ratios in air entering the stratosphere, with a 
peak-to-peak amplitude of ~3.2 ppmv. 

Both the seasonal cycle and interannual variations in the 
growth rate are potentially large enough to cause substantial 
error in estimates of mean age from /Lag> as noted in a model- 
ing study by Hall and Prather [1993]. For example, air with a 
true mean age of 0 year on July 1, 1992, would have a C0 2 
mixing ratio of 357.35 ppmv, according to our continuous rep- 
resentation of the C0 2 stratospheric boundary condition. How- 
ever, (2) would give a value of -1.2 year for /lag- Similarly, 
/lag for air with a true mean age of 0 year on July 1, 1993, 
(C0 2 = 357.85 ppmv) would be -0.6 year. If the C0 2 growth 
rate were constant at 1.47 ppmv yr , /Lag would equal 
-1.1 year for air with a true mean age of 0 year on July 1 of any 
year, since that is the date of the seasonal maximum in C0 2 
mixing ratios and the peak-to-peak amplitude of the seasonal 
cycle is ~3.2 ppmv. 

As air moves through the stratosphere, air parcels with dif- 
fering transit times since crossing the tropical tropopause mix, 
and the effects of seasonal and interannual variations in C0 2 
on mean age are attenuated. The amount of attenuation de- 
pends on the width of the age spectrum [Hall and Plumb , 
1994; Andrews et ai> 1999, 2001 ]. The examples above for air 
with a true mean age of 0 year represent a special case, where 
the age spectrum is a delta function with a spectral width of 0 
year. 

We have shown elsewhere that for air younger than -3 years 
the comprehensive data set for C0 2 from 1992 to 1998 pro- 
vides sufficient information to derive detailed empirical age 
spectra [Andrews et ai , 1999, 2001]. The seasonal cycle and 
growth rate changes were used in those analyses to estimate 
the width and shape of the age spectrum. We can thus use the 
first moment of the derived age spectra to provide accurate 
estimates of / as a function of N 2 0 for air younger than 
3 years. As shown in section 4.1, mean ages calculated from 
the age spectra are in good agreement with C0 2 lag times av- 
eraged over the 1992-1998 data sets for latitudes with frequent 
coverage. 

Little is known, however, about the age spectral width for 
older air, and it is not immediately obvious that using a single 
average growth rate from 1980 to 1998 will provide accurate 
estimates of the mean age. The growth rate calculated over 
various averaging periods is shown in Figure lb. For example, 
in mid- 1995 the “instantaneous’' growth rate of the deseasonal- 
ized C0 2 boundary condition is 2.2 ppmv yr 1 , while the aver- 
age calculated over the preceding 5 year period (mid- 1991 to 
mid- 1995) is only 1.1 ppmvyr 1 . These values bracket the 
long-term average growth rate of 1 .47 ppmv yr 1 . For a 5 ppmv 
lag in C0 2 relative to the boundary condition, the mean age 
corresponding to these growth rates would be 2.3, 4.5, and 
3.4 years, respectively. Intuitively, an average growth rate 
computed over a shorter (longer) interval seems more appro- 
priate for determining mean ages for younger (older) air, but 
implementing such a strategy is complicated. Fortunately, our 
analysis shows that errors associated with the propagation of 



seasonal and interannual variations in tropospheric C0 2 into 
the stratosphere are small (<0.3 year) for air with T > ~2 years, 
suggesting that the stratospheric age spectra are sufficiently 
wide to attenuate C0 2 fluctuations around the long-term trend 
during the period over which the sampled air entered the 
stratosphere. 

Figure 2 illustrates another potential complication for deriv- 
ing stratospheric mean ages from measurements of C0 2 . C0 2 
mixing ratios in air extracted from ice cores [Etheridge et ai , 
1996] provide a measure of atmospheric C0 2 levels prior to 
1958, when the first routine measurements of atmospheric C0 2 
were begun [Keeling and Whorf ’ 1994]. These data indicate a 
dramatic change in the C0 2 growth rate between 1960 and 
1980 (Figure 2a). In order to investigate whether such a 
change in the slope of the C0 2 boundary condition could sig- 
nificantly affect the mean age calculation, theoretical age spec- 
tra with a mean age of 5 years and varying width (Figure 2b) 
were convolved with a boundary condition comprised of the 
long-term trend shown in Figure 2a plus a seasonally varying 
component (as in Figure 1). 

As shown in Figures 2c and 2d, age spectra with long expo- 
nentially decaying tails pick up enough very old air to produce 
significant deviations from the C0 2 time series that would 
result from a constant growth rate. However, after 1984 the 
calculated mean ages are within 0.5 year of the true mean age 
for all but the widest test spectrum, which has a modal time 
(time of maximum probability) of 0.3 year and a tail so long 
that 13% of the air has been in the stratosphere more than 10 
years (6% for more than 20 years). This anomalous age spec- 
trum would require rapid transport of a large amount of air 
from the troposphere to the measurement location. However, 
mean ages of 5 years were found only for midlatitude and high- 
latitude balloon flights above 25 km and at ER-2 altitudes in 
filaments of air that had descended in the polar vortices. 
Large-scale transport of air from the tropical tropopause to 
altitudes >25 km on timescales of a few months is not consis- 
tent with current understanding of stratospheric transport. 
Thus a change in slope of the C0 2 stratospheric boundary con- 
dition prior to 1980 is unlikely to cause large errors in mean 
age estimates from measurements made from 1992 to 1998. 
For spectra with A 2 fT < 5 years, errors in the calculated mean 
age are <10%. Note that before 1980, the change in the C0 2 
growth rate is evident in the calculated C0 2 time series even 
for relatively narrow age spectra. This should be considered 
when comparing recent measurements with earlier strato- 
spheric C0 2 data, as in section 4.3.2. Tests using age spectra 
with mean ages of 4 and 7 years yielded similar results. 


4. Results and Discussion 

4.1. Results From the ER-2 Aircraft Data 

Boering et al [1994, 1996] showed that, below 21 km, 
stratospheric C0 2 and NiO are compactly related over a wide 
range of latitudes. N 2 0 is photolyzed rapidly at high altitudes 
but has a relatively long lifetime (years) below ~30 km. Thus, 



on average, N 2 0 mixing ratios decrease with altitude through- 
out the stratosphere. There are significant motions of air asso- 
ciated with planetary scale waves that displace tracer isopleths 
vertically but do not affect relationships among long-lived trac- 
ers [Ehhalt et al., 1983]. To remove this source of apparent 
variance, we use N 2 0 as a convenient vertical coordinate for 
analysis of C0 2 data. The spatial distribution of mean age is 
then obtained using mean distributions of N 2 0 with potential 
temperature and equivalent latitude [Strahan et al., 1999]. 
Equivalent latitude is a coordinate based on potential vorticity, 
which is nearly conserved on timescales typical of reversible 
wave transport [e.g., Nash et al., 1996]. By analyzing the C0 2 
data using /Ia G > we can explore the consistency of the relation- 
ship between C0 2 and N 2 0 from deployment to deployment 
and from year to year, which is otherwise complicated by the 
long-term increase in C0 2 mixing ratios. 

Figure 3 shows r LAG plotted as a function of N 2 0 for all of 
the ER-2 data collected from 1992 to 1998. The data span 
latitudes from 70°S to the North Pole with data in both hemi- 
spheres for all seasons except Southern Hemisphere summer. 
N 2 0 data were corrected for tropospheric growth over this pe- 
riod as described above. The la uncertainty on the average 
r L ag is < 0.3 year for air with 70 ppbv < N 2 0 < 275 ppbv. 
The small standard deviation in r LA0 for a given value of N 2 0 
provides evidence that age spectra in the lower stratosphere are 
wide enough to dampen seasonal variations and fluctuations in 
the C0 2 growth rate that occurred over the time period that the 
sampled air entered the stratosphere. Air with N 2 0 = 
275 ppbv has a mean age of 2.0 years, corresponding to a C0 2 
difference of 2.94 ppmv relative to the boundary condition, 
assuming a growth rate of 1 .47 ppmv yr . The observed 
la variation of -0.3 year in the calculated mean age corre- 
sponds to a variance in the average C0 2 growth rate of 
+0.2 ppmv yr 1 . This represents an upper limit, neglecting 
uncertainty in the measured C0 2 and N 2 0 mixing ratios. For 
older air, with larger C0 2 differences and smaller uncertainty 
in Tlag (< 0.2 year), the C0 2 growth rate is even more tightly 
constrained. Air with /Lag = 5 years (N 2 0 = 1 10 ppbv) re- 
flects a nearly invariant average C0 2 growth rate of 
1.47 + 0.05 ppmv yr 1 . Uncertainty in r LAG increases for 
N 2 0 < 80 ppbv, but these values represent relatively few data 
points, which were obtained in filaments of polar vortex air 
that had mixed to varying degrees with air from midlatitudes. 

Figure 4a shows r LAG plotted as a function of N 2 0 for data 
collected from the ER-2 aircraft within 2.5° of the NASA 
Ames Research Center in Moffett Field, California (37.4°N) 
from 1992 to 1998 during each of the four seasons. Note that 
there is very little seasonal variation in this relationship for 
200 ppbv < N 2 0 < 275 ppbv. The corresponding standard 
deviations indicate that interannual variations are also small 
(Figure 4b). For air with N 2 0 > 275 ppbv (corresponding to 
altitudes < -19 km at this latitude), the influence of the C0 2 
seasonal cycle propagating from the troposphere via the tropi- 
cal stratosphere is evident: artificially young (old) mean ages 
are calculated in summer (winter) when the observed C0 2 con- 
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centration in the extratropical stratosphere is at its maximum 
(minimum). Andrews et aL [2001] showed that the CO 2 sea- 
sonal cycle is still discernible for air with N 2 0 = 215 ppbv but 
that the peak-to-peak amplitude is reduced to 0.8 ppmv, corre- 
sponding to small errors in the mean age (< ±0.3 year). The 
seasonal variations in /"observed for air with N 2 0 < 200 ppbv 
are apparently due not to the propagation of tropospheric sea- 
sonal variations, but rather to the mixing of remnants of the 
polar vortex with air from lower latitudes during spring 
[Waugh et al. , 1997]. 

When averaged over all seasons and years, the relationship 
between /Lag and N 2 0 is remarkably invariant over a large 
range of latitudes (Figure 4c). Nevertheless, estimates of the 
mean age from /Lag for a particular flight are subject to sig- 
nificant error for air with N 2 0 > 275 ppbv due to seasonal and 
interannual variations in C0 2 mixing ratios propagating from 
the tropical tropopause, as evidenced by the sharp increase in 
the standard deviation of the average /Lag at all latitudes (Fig- 
ure 4d). For the individual deployments, typical (Northern 
Hemisphere) NH midlatitude C0 2 standard deviations in a 
10 ppbv N 2 0 bin centered at 295 ppbv were < 0.4 ppmv [An- 
drews et al , 2001], corresponding to uncertainties in /Lag of 
< 0.3 year. Thus the increased standard deviations seen in 
Figures 4b and 4d are not explained by increased measurement 
uncertainty at high values of N 2 0. 

A method for estimating /"for a particular flight that is not 
sensitive to errors associated with seasonal and interannual 
variations in C0 2 mixing ratios is highly desirable. Andrews 
et al [2001] showed that for NH midlatitudes, empirical age 
spectra can be derived from the time series of C0 2 observations 
for air with N 2 0 > 235 ppbv, where seasonal and interannual 
oscillations in C0 2 propagating from the troposphere provide 
information about transport processes occurring on multiple 
timescales. They found that the data are consistent with a sea- 
sonally invariant relationship between mean age and N 2 0. 
This relationship is essentially equivalent to the cubic fit to the 
averaged /Lag:N 2 0 correlation shown in Figure 3, 
within ± 0.05 year. 

Thus for observations at NH midlatitudes we can apply the 
relationship between mean age from the empirical age spectra 
r SPE c and N 2 0 to estimate mean ages from measured N 2 0 for 
air with N 2 0 > 275 ppbv, thereby avoiding complications 
caused by the seasonal cycle and variability in the C0 2 growth 
rate. For N 2 0 data with uncertainty < 2.5%, this method pro- 
vides mean ages with an uncertainty of < 0.5 year. Figure 4c 
suggests that the average /spec'N 2 0 relationship provides 
mean ages with errors of < 0.5 year regardless of latitude at 
ER-2 altitudes. (Whether this result extends to higher altitudes 
is considered in section 4.2.) 

To further explore whether the relationship between /spec 
and N 2 0 derived from midlatitude data applies in the tropics, 
we compare to mean ages from tropical age spectra derived by 
Andrews et al [1999] (Figure 5). Mean ages from the tropical 
age spectra agree to within 0.3 year with those predicted by 
applying the midlatitude /"spec^O relationship to tropical 



N 2 0 data, except in NH late winter/early spring, when the 
mean ages from the tropical age spectrum analysis are young- 
est. The observed tropical N 2 0 mixing ratios were slightly 
lower in March 1994 and February 1996 than at other times of 
year, which should correspond to older air according to the 
relationship between /spec an d N 2 0 derived for midlatitudes. 
As noted by Andrews et al. [1999], the seasonality in the de- 
rived tropical age spectra relies heavily on the results for a 
single profile (February 1996). However, the variations in the 
tropical N 2 0 mixing ratios are of the order of the accuracy of 
the ATLAS instrument, so it is not clear whether the apparent 
differences in the seasonal mean ages are significant. In any 
case, the maximum error incurred by applying the midlatitude 
r SPE c:N 2 0 relationship in the tropics is 0.5 year. Although this 
amounts to a difference of nearly 100 % for these young ages, 
knowledge of the tropical mean age to better than 0.5 year 
places a tight constraint on models of the stratospheric circula- 
tion. 

We showed above that r L AG provides a precise estimate of 
the mean age for air with N 2 0 < 275 ppbv and that we can use 
the relationship between /"spec and N 2 0 derived for midlati- 
tudes to estimate mean age for air with N 2 0 > 275 ppbv, re- 
gardless of latitude. While this approach produces the correct 
average distribution of mean age, it creates a small discontinu- 
ity in the mean age at N 2 0 = 275 ppbv for most flights. To 
avoid this, we introduce the following functional form for r 
with a smooth transition from /spec 1 ° ^ AG ‘ 

r= /Lag = t - ( 1749.88 + 0.6824 1 (C0 2 - ( 1 .725 - CH 4 ))) 
N 2 0 < 235 ppbv 


r= 0.025(275 - N 2 0)/L ag + 0 - 0.025(275 - N 2 0 )) /spec 
235 ppbv < N 2 0 < 275 ppbv 


r= r s pec =0.0566(313 - N 2 0 ) - 0.000195(313 - N 2 0) 2 
N 2 0 > 275 ppbv. (3) 

Here the units of C0 2 and CR, are ppmv, and N 2 0 has units of 
ppbv. Equation (3) produces a smooth profile of mean age for 
each flight. Errors associated with the persistence of seasonal 
and interannual variations in the growth rate of C0 2 are < 0.5 
year. In the following discussion the symbol F refers to mean 
ages calculated using (3). 

Figure 6 a shows the latitudinal distribution of r at 
20 ± 0.5 km, corresponding to the cruise altitude of the ER-2 
aircraft. The shaded symbols represent 10-s averages of the 
data along the flight track, and the large solid symbols with 
error bars represent the median value in overlapping 5° latitude 
bins. Note the large meridional gradients in the subtropics of 
both hemispheres. Similar gradients have been observed for 
other tracers, for example, NQ/O3 [Murphy et al . , 1993]. This 
feature is a sensitive indicator of meridional exchange between 
the tropics and higher latitudes [Plumb, 1996] that most mod- 
els of the stratosphere currently fail to simulate [Park et al . , 
1999; Hall et al, 1999]. 



The latitudinal variation of r is compared in Figure 6b with 
that calculated by simply averaging /Lag- Since NH data are 
available for many seasons over a period of several years, r 
and /Lag produce nearly identical results, as expected based on 
Figures 3 and 4. There is good agreement between /"and r , a g 
for the middle and high latitudes of both hemispheres, which 
implies that the corresponding age spectra are sufficiently 
broad to dampen the effects of seasonal and interannual varia- 
tions in C0 2 during the time period over which the data were 
acquired. 

In the southern tropics and subtropics, however, /Lag is as 
much as 0.8 year older than F The reason for this difference is 
twofold. First, we do not have Southern Hemisphere (SH) data 
for December - February (Figure 6c), the time of seasonal 
maximum in C0 2 concentrations in the tropics at 20 km. The 
seasonal maximum in C0 2 corresponds to artificially young 
ages, and since these data are lacking, the average is skewed 
toward older ages. Second, as discussed above, interannual 
variations in the C0 2 growth rate propagating into the strato- 
sphere can cause significant (for air with N 2 0 > 275 ppbv) 
bias in /Lag, which is calculated assuming a linear growth 
rate. The SH data shown in Figure 6a were acquired from 
March to October 1994; thus the data from the equator to 20°S, 
with values for /"of 1 ± 0.5 year, represent air that entered the 
stratosphere between early 1993 and early 1994. According to 
Figure 1, the annual mean (i.e., deseasonalized) C0 2 mixing 
ratios in air entering the stratosphere were lower than the linear 
fit used to calculate /Lag (equation (2)) by 0.25 to 1 ppmv 
throughout this time period, which corresponds to a bias to- 
ward older ages of 0.1 - 0.7 year and is therefore consistent 
with the size of the observed discrepancy. Based on the NH 
data, /"and /Lag should converge if more seasons and years are 
sampled in the SH stratosphere. 

The only significant seasonal variation in the average r at 
20 km is found in the NH subtropics where mean ages during 
December-February are older by -1 year than for March- 
August and older by ~2 years than for September-November 
(Figure 6c). There is also some variation in the SH subtropics, 
but the data for this region are temporally too sparse to inter- 
pret, as discussed above. Subtropical variations in mean age 
may be associated with the seasonal migration of the boundary 
of the tropical upwelling region. Using observations of the 
distribution of stratospheric aerosols, Grant et al. [19961 found 
that the northern boundary of the tropics did not extend north 
of 20°N from November to March 1991, while there were fre- 
quent excursions farther north during the rest of the year. A 
mean age maximum during NH winter is consistent with this 
migration. Note that we would also expect to see significant 
seasonal variations at high latitudes associated with the polar 
vortices, but this analysis does not include data for the south 
polar region during southern summer or for the north polar 
regions during northern winter. 

In order to obtain further information about the spatial dis- 
tribution of mean age, we used the polynomial relationship 
between average /Lag an d N 2 0 shown in Figure 3 to generate 



distributions of mean age as a function of equivalent latitude 
and potential temperature, by mapping this relationship onto 
the seasonally resolved N 2 0 climatology of Strahan et ai 
[1999] (Figure 7). This function gives values of T that are 
within 10% of the value calculated using (3) for all but 2% of 
the individual measurements and that are within 15% of the 
value calculated using (3) for all but 0.2% of the data. The 
resulting distributions provide a comprehensive picture of 
mean age in the lower stratosphere that can easily be compared 
with output from numerical models. Contours of mean age are 
relatively symmetric about the equator during NH spring, and 
old air descending in the polar vortex can be seen during win- 
ters in both hemispheres. Note that our analysis did not in- 
clude any data north of ~60°N for the December-February pe- 
riod, but high-latitude winter data were included in the N 2 0 
climatology from earlier deployments. Thus we have assumed 
that the T?N 2 0 relationship is valid in the polar vortex at these 
altitudes. Data obtained in the 1999/2000 Arctic vortex during 
the Sage III Ozone Loss and Validation Experiment (SOLVE) 
indicates that this is valid. 

4.2. Extension to Higher Altitudes 

The OMS flights produced a unique set of in situ observa- 
tions of C0 2 , N 2 0, CH 4 , O3, H 2 0, SF 6 , CFC-1 1, CFC-12, and 
Halon-121 1 up to 32 km. Data were obtained in the tropics 
and at NH middle and high latitudes from 1996 to 1998 during 
several seasons. In this section, vertical profiles of mean age 
are presented for each of the balloon flights and the relation- 
ships between mean age and N 2 0 and CFC-1 1 above ER-2 
altitudes are explored. 

We first consider whether (3) can be used to obtain mean 
ages for the balloon profiles. This should not present a prob- 
lem at middle and high latitudes, since air with N 2 0 
> 275 ppbv was not encountered above 21 km, the maximum 
altitude of the ER-2. However, in the tropics, air with 
N 2 0 > 275 ppbv was sampled up to 23 km. Figures 8 a and 8 b 
show /Lag as a function of N 2 0 for each extratropical and 
tropical balloon flight. Average relationships between /lag 
and N 2 0 for the tropics and extratropics were produced from 
the OMS data by calculating the average relationship for each 
flight, then averaging the results according to latitude (Figure 
8 c). The February and November results were weighted 
equally in determining the average tropical relationship. The 
polynomial relationship derived from ER-2 data is shown for 
comparison in Figures 8 a and 8 c. The "reference correlation 
of Strunk et ai [2000], derived from C0 2 and N 2 0 data in 
whole air samples from 10 middle and high-latitude balloon 
flights, is also shown in Figure 8 a; agreement is relatively 
good, although we obtain ages that are 10 - 20 % older than their 
reference correlation for air with 275 ppbv > N 2 0 > 175 ppbv. 
Part of the discrepancy may be due to the fact that Strunk et al. 
calculate mean age relative to the surface rather than the tropi- 
cal tropopause, but this should amount to a difference of only 2 
months and would appear as an offset over the entire range of 
N 2 0 values. 



For air with N 2 0 > 275 ppbv, the agreement between the 
average relationship from the ER-2 data and the average /" ^ AG 
for both the tropical and the extratropical balloon flights is 
somewhat surprising, given the limited number of balloon 
profiles. Errors in the averaged r LAG resulting from seasonal 
and interannual variations in C0 2 are fortuitously negligible 
for the tropical balloon flights due to fortunate timing. For 
example, the age spectrum analysis of the tropical ER-2 data 
[Andrews et ai , 1999] shows that the vertically propagating 
C0 2 seasonal cycle in February is roughly 180° out of phase 
with respect to November below 460 K, so we expect that er- 
rors in /" LAG associated with the persistence of the seasonal 
cycle will largely cancel when data from these months are av- 
eraged. At higher potential temperatures, error due to the sea- 
sonal cycle would be < 0.4 year, even if there were no cancella- 
tion, because the peak-to-peak amplitude is small [Andrews et 
al y 1999]. Also from 1995 to 1998 the long-term trend in C0 2 
mixing ratios entering the stratosphere is nearly indistinguish- 
able from the linear fit used to calculate /I AG . Thus for a 
“young” air parcel that does not contain a large amount of air 
that entered the stratosphere prior to 1995, errors due to inter- 
annual variation in the C0 2 growth rate are small for all of the 
OMS flights. Older air parcels are likely to have wide enough 
age spectra to dampen errors due to fluctuations in the C0 2 
growth rate. We therefore conclude that it is legitimate to use 
(3) to generate smooth profiles of mean age for each of the 
OMS balloon flights to date. Note, however, that uncertainty 
in /"depends on the accuracy of available N 2 0 data. Uncer- 
tainty in individual values of /Tor OMS flights can be as large 
as ± 1 year because available N 2 0 data typically has an esti- 
mated accuracy of 5%. 

The ER-2 and OMS data demonstrate that the relationship 
between /"and N 2 0 is remarkably invariant over a wide range 
of latitudes and altitudes. In contrast, there is significant lati- 
tudinal variation in the relationship between r and CFC-11 
(Figure 8d). The existence of a compact relationship between 
two species with different chemical lifetimes provides infor- 
mation about transport timescales [Plumb and Ko , 1992]. 
Evidently, at ER-2 and OMS altitudes, horizontal exchange 
across the subtropics occurs rapidly compared to the lifetime of 
N 2 0 (a few years at 30 km), but photochemical loss is faster 
than meridional transport for CFC-1 1, which has a lifetime of 
a few months at these altitudes. We would expect the invari- 
ance of the relationship between F and N 2 0 to break down 
where the lifetime of N 2 0 becomes very short, for example, at 
higher altitudes in the tropics. The ability to simulate relation- 
ships among tracers with varying lifetimes is a key test of 
transport parameterizations in models of the stratosphere. 

Mean ages calculated using (3) are shown versus potential 
temperature in Plate 1 for the OMS flights. Tropical ages are 
~2 years younger than extratropical ages on the same isentrope 
over most of the potential temperature range. Neu and Plumb 
[1999] showed that in a simple model this difference is a fun- 
damental property of the circulation that depends primarily on 
the vertical velocity in the tropics. These results represent the 



first observations of the mean age difference between the inner 
tropics and midlatitudes above -20 km. Evidence of 
large-scale mixing events can also be seen in the mean age 
profiles in Plate 1. As discussed below, remnants of the 
1996/1997 Arctic polar vortex are evident in the June 1997 
profile from Fairbanks, Alaska (65 °N) near 520 and 625 K. 
Filaments of extratropical air appear in both of the Novem- 
ber 1997 profiles from Juazeiro do Norte, Brazil (7°S), indicat- 
ing that air at that latitude was not completely isolated from 
midlatitudes at that time. 

The May 1998 OMS flight from Fort Sumner, New Mexico 
(35°N) is characterized by a region of nearly constant mean 
age above -550 K (-22.5 km). SF 6 data from the September 
1996 flight indicate a similar feature above -650 K (-26 km). 
Unfortunately, there are no C0 2 data available above 640 K 
(-25.5 km) for that flight. A region of nearly constant mean 
age is present in almost every midlatitude balloon profile of 
C0 2 and SF 6 to date, starting between 20 and 26 km and ex- 
tending to the maximum altitude of the profile (-32 km), with 
typical mean ages of 4-6 years in this region [Bischof et ai, 
1980, 1985; Schmidt and Khedim , 1991; Nakazawa et al, 
1995; Harnisch et al. 1998]. We present a detailed compari- 
son of these measurements and consider implications for mod- 
els of transport in section 4.3.2. 

Although the relationship between r and N 2 0 has been 
shown to be remarkably compact, careful examination of data 
from the nearly coincident balloon and ER-2 flights of June 30, 
1997, reveals subtle but significant differences in the N 2 0:/" 
relationship at ER-2 altitudes versus higher altitudes (Plate 2). 
Using data from ATMOS, Michelsen et al. [1998] suggested 
that distinct relationships between mixing ratios of N 2 0 and 
CH 4 can be used to identify air with tropical, midlatitude, or 
polar vortex character. Herman et al. [1998] showed that the 
N 2 0 :CH 4 relationships observed during these OMS flights are 
in good agreement with the ATMOS relationships and that air 
with both midlatitude and vortex character was sampled during 
the June 30, 1997, balloon flight. Air characterized by the 
'Vortex-like" N 2 0 :CH 4 relationship corresponds to distinct 
features in the tracer profiles, the sharp minima in C0 2 
(maxima in mean age) at -520 K and 640 K (Plate 2a). The 
presence of such sharp features was a surprise given that the 
vortex had broken up several months earlier. Note that the 
ER-2 only sampled air with low N 2 0 values in the vortex rem- 
nants, while the balloon sampled air with low N 2 0 both in the 
vortex remnants and aloft. From 10 km to the maximum alti- 
tude of the ER-2, the median difference between the ER-2 and 
balloon C0 2 data was 0.04 ppmv, well within the combined 
uncertainties of the instruments (Daube et al., submitted 
manuscript, 2001 ). 

The relationship between C0 2 and N 2 0 for this flight is 
shown in Plate 2b for air with r > 3.5 years (C0 2 > 
358 ppmv). Balloon data classified as "midlatitude-like" ex- 
hibit a slightly sigmoidal (S-shaped) correlation, while the 
correlations corresponding to the vortex remnants are signifi- 
cantly straightened. One interpretation of the data in Plate 2b 



is that the correlations for the vortex remnants reflect downwel- 
ling of midlatitude-like air from high altitudes in the polar 
vortex followed by mixing with midlatitude air at lower alti- 
tudes, either as midlatitude air was entrained into the vortex 
[Plumb et al. , 2000] or after the vortex had broken up [e.g., 
Waugh, 1997], resulting in a straightening of the relationship 
across the "concave down" portion of the high-altitude correla- 
tion. However, data from the April 26, 1997, ER-2 flight are 
included in Plate 2b to show that the correlations in the vortex 
remnants in July could also be interpreted as a straightening of 
the relationships across the "concave up" correlation observed 
near the edge of the polar vortex in April. For comparison, 
Plate 2c shows the relationship between r and N 2 0 for the 
same data. Because mean age is not conserved in the strato- 
sphere (having a uniform source of 1 yr/yr), in general it is 
more appropriate to examine particular mixing events using 
C0 2 corrected for CH 4 oxidation, which is conserved. Ac- 
counting for the source of mean age would be equivalent to 
adding 0.18 year to the mean ages for April 1997 in Plate 2c. 

43. Comparison With Other Estimates of Mean Age 

43.1. Mean ages from simultaneous measurements of 

SF 6 . We first consider mean ages derived from measurements 
of SF 6 obtained concurrently with the C0 2 data used in this 
analysis. Figure 9 shows a scatterplot of all the mean ages 
derived from C0 2 and SF 6 measurements obtained using the 
ER-2 aircraft (-3000 points). Mean age estimates are derived 
from stratospheric SF<; mixing ratios using the expressions of 
Volk et al. [1997]. The boundary condition is the global mean 
surface time series of SF 6 [Geller et al., 1997] with an imposed 
delay of 0.8 year representing the transit time for surface air to 
reach the tropical tropopause. The delay time is chosen to give 
the best agreement with measured tropical tropopause concen- 
trations where available. 

There is excellent agreement between the C0 2 and SF 6 es- 
timates for mean ages up to -3 years, while for older air the 
SF 6 ages are - 20% older. Results are similar for OMS data, 
with mean ages derived from SF 6 10-20% older than the corre- 
sponding mean ages from C0 2 [see Park et al., 1999; Hall et 
al., 1999]. Strunk et al. [2000] reported somewhat better 
agreement (± 0.5 year, - 10 %) between mean ages derived from 
SF(, and C0 2 in whole air samples obtained during balloon 
flights at midlatitudes and high latitudes referenced to the sur- 
face rather than the tropical tropopause. 

There are several possible explanations for differences be- 
tween mean ages derived from SF 6 and C0 2 . The growth rate 
for SF 6 is slightly nonlinear, and a correction factor based on 
the width of age spectra in the lower stratosphere from models 
was applied to calculate the SF^ ages [Volk et al., 1997]. Errors 
in the correction could lead to a systematic bias that would be 
largest for old ages. Note that Volk et al. [1997] claimed only 
that the correction was valid for the lower stratosphere 
(< 21 km), although here and in the work of Strunk et al. 
[ 2000 ] the correction was applied to balloon data obtained up 
to -30 km. Also uncertainty in the stratospheric boundary 



condition for SF 6 gives rise to uncertainty of the order of 
± 0.5 year in derived mean ages. Finally, there is the possibil- 
ity that mesospheric loss of SF 6 leads to artificially old mean 
ages from SF 6 observations [Hall and Waugh , 1998]. Some 
evidence for the latter theory is found in the data from the 
June 30, 1997, OMS flight from Fairbanks, Alaska (65°N). 
Mean ages from C0 2 were oldest at the top of the profile, while 
mean ages from SF 6 were significantly older in remnants of the 
polar vortex [see Hall et ai, 1999, Figure 5], consistent with 
the idea that SF 6 is destroyed at high altitudes. A similar 
anomalous discrepancy between mean ages from SF 6 and C0 2 
in whole air sampled during a balloon flight over southern 
France on June 23, 1997, has been attributed to mesospheric 
loss of SF 6 [Strunk et ai , 2000]. Preliminary analysis of data 
from the 1999/2000 Arctic vortex during the SOLVE cam- 
paign also indicates substantial mesospheric loss of SF 6 (F. L. 
Moore, manuscript in preparation, 2001). 

4,3.2. Comparison with previous measurements of C0 2 
in the middle stratosphere. Schmidt and Khedim [1991] 
compiled a list of C0 2 measurements from balloon-borne cryo- 
samplers through 1990 that includes 21 flights at middle or 
high latitudes in the Northern Hemisphere. All of these pro- 
files were characterized by a region of nearly constant C0 2 
mixing ratio at the top of the profile. They reported an average 
mean age corresponding to this region of 5.6 ± 1.1 years, in 
good agreement with the mean ages from the OMS flights over 
Fort Sumner, New Mexico (35°N) and Fairbanks, Alaska 
(65°N), although with relatively large uncertainty. A time 
series of the average mixing ratio corresponding to the region 
of near-constant C0 2 for each of the profiles given by Schmidt 
and Khedim [1991] is shown in Figure 10 along with C0 2 data 
from the same region for OMS flights presented here. Data 
from Nakazawa et ai [1995] and Harnisch et ai [1998] are 
also shown. 

The data presented in Figure 10 have not been corrected for 
stratospheric C0 2 production by CH* oxidation. During OMS, 
CFLjin the midlatitude middle stratosphere was -1000 ppbv on 
average, and the tropospheric value was -1725 ppbv. Thus the 
correction applied to the C0 2 data would be 1.725 - 
1.000 = 0.725 ppmv, or 42% of the tropospheric CH* mixing 
ratio, and neglecting to account for stratospheric C0 2 produc- 
tion by CH4 oxidation would result in mean ages that are too 
young by approximately (0.725 ppmv)/(1.47 ppmvyr 1 ) 
= 0.5 year. In 1975 the tropospheric CH4 mixing ratio was 
-1440 ppbv, implying a correction to the midlatitude middle 
stratospheric C0 2 data of 0.60 ppmv, according to the relation- 
ship between midstratospheric and tropospheric CH 4 observed 
during OMS. Differences in the C0 2 correction due to the 
tropospheric trend in CH4 over this time period are therefore 
likely to be < 0.2 ppmv, unless the rate of stratospheric CH4 
oxidation has changed substantially since 1975. 

When corrected for CH4 oxidation, the C0 2 data from the 
extratropical OMS flights correspond to a lag time with respect 
to the stratospheric boundary condition of 5.0 ± 0.5 years. 
This value is notably consistent with the Nakazawa et ai 
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[1995] data obtained from 1985 to 1991 and with the Bischof 
et al [1985] and Volz et al [1981] data obtained from 1976 
tol984 (as reported by Schmidt and Khedim [1991]), although 
in general the earlier data have larger uncertanties. Data for 
the OMS flight of September 21, 1996, are not included since 
there are no C0 2 data above ~25 km for that flight, but analy- 
sis of concurrent observations of SF 6 indicates that the mid- 
stratospheric mean age is similar to that obtained for the May 
18, 1998, flight. For the OMS and Nakazawa et al [1995] 
data, the error bars represent the standard deviation of the C0 2 
mixing ratio in the near constant region. For OMS the meas- 
urement precision is 0.1 ppmv, corresponding to the size of the 
symbols in Figure 10b. Error bars for Volz et a. I [1981] and 
Bischof et al. [1985] are as given by Schmidt and Khedim 
[1991]. 

The data reported by Schmidt and Khedim [1991] and by 
Harnisch et al. [1998] have considerably more scatter than the 
other data sets. The error bars in Figure 10b on the Schmidt 
and Khedim data are as given in that study, but a later paper by 
that group [Strunk et al. , 2000] stated that the precision on the 
earlier measurements was ±1 ppmv. Schmidt and Khedim 
attributed the variability in their calculated ages to real changes 
in the large-scale stratospheric circulation, but that interpreta- 
tion seems inconsistent with the observations of Nakazawa et 
al. [1995] during the same time period, which show little vari- 
ability. 

The C0 2 values shown for Harnisch et al. [1998] were ob- 
tained by averaging their midlatitude data for altitudes greater 
than 20 km, and the error bars represent 1 standard deviation 
of the mean. Their stated measurement precision is ± 
0.2 ppmv. Interestingly, they report mean ages from SF 6 for 
both of these flights that are consistent with the 5.0 ± 0.5 year 
mean ages implied by the ensemble of C0 2 data. Moreover, 
while their SF 6 profiles are characterized by a smooth decrease 
with a region of near-constant mixing ratio at the top of the 
profile, like the C0 2 and SF ft profiles discussed above, their 
C0 2 profiles are noisy. They also reported N z O mixing ratios, 
which facilitates comparison of their observations from 1993 
and 1995 at 20 km (215 ± 5 pbbv N 2 0 for both flights) with 
C0 2 data from the ER-2 from 1992 to 1998 and from the 
OMS balloon flights (Figure 11). The height of the symbols 
reflects the la standard deviation in the average of the ER-2 
and OMS C0 2 data. The ER-2 and OMS data show a linear 
increase with time on the 215 ppbv N 2 0 isopleth, consistent 
with the average tropospheric trend, but Harnisch et al [1998] 
observed a decrease in corresponding C0 2 mixing ratios be- 
tween 1993 and 1995. This decrease is even larger at other 
altitudes (1-2 ppmv), despite very similar N 2 0 profiles for the 
two flights. It is difficult to imagine a plausible transport sce- 
nario capable of explaining these discrepancies. Strunk et al 
[2000] reported good agreement between mean ages calculated 
from SF 6 and C0 2 in whole air samples from the Arctic and 
midlatitude stratosphere and argued that the C0 2 mixing ratios 
reported by Harnisch et al [1998] may have been affected by 
contamination or sample degradation. 


The consistency between our results and those of Nakazawa 
et al [1995], Bischof et al [1980, 1985], and Volz et al 
[1981] suggests that the mean residence time of air in the mid- 
dle stratosphere has remained within 1 year of its current value 
over the past 25 years. However, as noted above, we would 
expect C0 2 mixing ratios to have been influenced by the dra- 
matic change in the C0 2 growth rate that occurred between 
1960 and 1980. Figure 2d indicates that air with a mean age 
of 5 years would have C0 2 mixing ratios at least 1.5 ppmv 
higher than predicted for a constant C0 2 growth rate, if the real 
age spectrum resembles any of the test spectra. The fact that 
C0 2 mixing ratios were not measurably elevated in the early 
flights relative to the long-term trend suggests that the age 
spectrum for this region is rather narrow, although it is possi- 
ble that transport rates have changed over this period in such a 
way as to mask a change in the slope of the C0 2 time series. 
The ability of models to reproduce the midlatitude midstrato- 
spheric C0 2 time series will provide a useful diagnostic. This 
analysis illustrates the potential value of a time series of high- 
accuracy C0 2 measurements as an indicator of changes in the 
circulation. In addition, future changes in the C0 2 growth rate 
may allow more detailed information about the age spectrum to 
be inferred from observations. 

The repeated observation of a region characterized by 
near-constant mean age in the midlatitude stratosphere also 
has implications for understanding the general circulation of 
the stratosphere. Bischof et al. [1985] found that profiles with 
this shape could not be reproduced using a simple time- 
dependent one-dimensional diffusion model and suggested that 
the shape may result from tropical upwelling followed by 
poleward transport in the middle and upper tropical strato- 
sphere. Nakazawa et al. [1995] noted that the profiles are 
consistent with more rapid poleward motion near the tro- 
popause than at higher altitudes, a transport scenario that has 
also been invoked to explain numerous observations from sat- 
ellites, balloons, and aircraft over the last decade. 

The shape of the midlatitude C0 2 profiles, as well as the 
observed mean ages in the midlatitude middle stratosphere, 
provide corroborative evidence for the validity of the bi modal 
age spectra derived for the midlatitude lower stratosphere from 
this time series of C0 2 observations by Andrews et al. [2001]. 
The younger peak in their derived age spectra (corresponding 
to air younger than 1 year) was thought to represent air trans- 
ported quasihorizontally from the tropical stratosphere, while 
the older peak (consisting of air that has been in the strato- 
sphere for 4 to 6 years) represents air descended from higher 
altitudes. The mean ages observed from 25 to 32 km during 
the OMS flights at middle and high latitudes of 5 ± 0.5 years 
are consistent with that region being the source of the air com- 
prising the older peak in the bimodal age spectra derived for 
lower altitudes. The separation of the two peaks in the age 
spectra requires a significant altitude region above 20 km 
where meridional transport is slow relative to higher and lower 
altitudes, consistent with the idea of a subtropical “barrier 1 ’ to 
horizontal mixing [e.g., Plumb , 1996]. 
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5. Conclusions 

We have derived accurate mean ages for stratospheric air 
from a comprehensive set of in situ observations of C0 2 , CH4, 
and N 2 0 obtained from 1992 to 1998 from aircraft and balloon 
flights. Detailed analysis of the uncertainty in the derived 
mean ages was possible because of the depth of the underlying 
dataset, characterized by extensive latitudinal coverage and 
frequent sampling. Errors associated with the seasonal cycle 
in C0 2 mixing ratios and interannual variation in the C0 2 
growth rate are < 0.5 year throughout the stratosphere (for r 
defined by (3) or for averaged 71 a g values) and <0.3 year for 
air with N 2 0 < 275 ppbv (for any r for a given C0 2 :N 2 0 
measurement). We conclude that stratospheric age spectra are 
broad enough to attenuate the effects of the seasonal cycle and 
fluctuations in the C0 2 growth rate over the time period during 
which the sampled air entered the stratosphere. 

A remarkably compact relationship exists between meas- 
ured N 2 0 and mean ages derived from C0 2 that is independent 
of latitude. In contrast, the relationship between CFC-11 and 
mean age exhibits substantial latitudinal variation. These rela- 
tionships give insight into timescales for exchange of air be- 
tween the tropics and higher latitudes. The robust relationship 
between mean age and N 2 0 allows a climatology of mean age 
to be derived from even more extensive observations of N 2 0 
during this period. 

The climatology of mean age can be used to evaluate trans- 
port in models used to assess the impact of anthropogenic per- 
turbations on stratospheric ozone concentrations. The distribu- 
tion of mean age with latitude and altitude provides detailed 
quantitative information about the general circulation of the 
stratosphere. At 20 km, sharp meridional gradients in the 
mean age are observed in the subtropics, further confirming 
earlier observations that exchange between the tropics and 
higher latitudes is significantly inhibited [e.g., Randel et al . , 
1993; Murphy et ai , 1993; Volk et al . , 1996; Plumb , 1996]. 
Between 20 and 30 km the average difference in mean age 
between the tropics and midlatitudes is -2 years, with slightly 
smaller differences at higher and lower altitudes. 

Vertical profiles of C0 2 mixing ratios at midlatitudes are 
characterized by a rapid decrease from the tropopause to a re- 
gion of near-constant mixing ratio in the middle stratosphere 
(-25-32 km), corresponding to mean ages of 5± 0.5 years, 
notably similar to the mean age of the older peak in the bi- 
modal age spectra derived for the midlatitude lower strato- 
sphere of Andrews et al. [2001]. This shape is consistent with 
the view that the meridional exchange is most rapid near the 
tropopause, with limited exchange from 20 to 30 km. Com- 
parison of C0 2 measurements obtained in the midlatitude mid- 
dle stratosphere since the mid-1970s with our balloon data 
shows that the mean age of air in this region has been rela- 
tively constant over the last 25 years, suggesting that continued 
monitoring of C0 2 at high altitudes could provide a baseline 
against which to measure changes in the stratospheric circula- 
tion. 
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Table 1. ER-2 Field Deployments 


Field Campaign 



Deployment Site(s) 

Dates 

Latitude Range 

SPADE 

Moffett Field, California 

Nov. 9-21, 1992 

20°N - 40°N 


Moffett Field, California 

April 19 to May 18. 1993 

14°N - 60°N 


Moffett Field, California 

Oct. 20- 25, 1993 

14°N - 60°N 

ASHOE/MAESA 

Moffett Field, California 

Feb. 2- 19, 1994 

35°N - 61°N 


Moffett Field, California 
Barber’s Point, Hawaii 
Nadi, Fiji 

Christchurch, New Zealand 

March 18 to April 15, 1994 

68°S - 38°N 


Christchurch, New Zealand 

May 24 to June 8, 1994 

69°S - 20°S 


Christchurch, New Zealand 

July 25 to Aug. 10, 1994 

69°S - 19°S 


Christchurch, New Zealand 

Oct. 3 to Nov. 4, 1994 

70°S - 60°N 

STRAT 

Moffett Field, California 

May 1 - 18, 1995 

14°N - 62°N 


Moffett Field, California 
Barber’s Point, Hawaii 

Oct. 18 to Nov. 9, 1995 

2°S - 59°N 


Moffett Field, California 
Barber's Point, Hawaii 

Jan. 23 to Feb. 15, 1996 

2°S - 53°N 


Moffett Field, California 
Barber's Point, Hawaii 

July 16 to Aug. 10, 1996 

2°S-61°N 


Moffett Field, California 

Sept. 12-21, 1996 

14°N - 61°N 


Moffett Field, California 
Barber’s Point, Hawaii 

Dec. 2- 19, 1996 

2°S - 61°N 

POLARIS 

Moffett Field, California 
Fairbanks, Alaska 

April 17 to May 13, 1997 

14°N - 90°N 


Fairbanks, Alaska 

June 22 to July 10, 1997 

45°N - 90°N 


Fairbanks, Alaska 
Barber’s Point, Hawaii 

Sept. 3 -25, 1997 

2°S - 90°N 

Flights-of- 

opportunity 

Edwards, California 

Nov. 10- 15, 1998 

1 1°N - 37°N 


Table 2. Observations of the Middle Stratosphere (QMS) 


Deployment Site 

Dates 

Latitude 

Altitude Range for CO 2 (km) 

Fort Sumner, New Mexico 

Sept. 21, 1996 

35°N 

9-26 

Juazeiro do Norte, Brazil 

Feb. 14, 1997 

7°S 

9-28 

Fairbanks, Alaska 

June 30, 1997 

65 °N 

9-32 

Juazeiro do Norte, Brazil 

Nov. 11 and 20, 1997 

7°S 

9-32 

Fort Sumner, New Mexico 

May 18, 1998 

35°N 

9-32 




Figure 1 . (a) Stratospheric boundary condition for CO 2 : symbols represent observations of 
CO 2 in air entering the stratosphere; the light curve is a continuous representation assuming a 
constant seasonal cycle and that the long-term trend can be represented using the dashed curve, 
which is the 12-month running mean of the average of surface station data from Mauna Loa, 
Hawaii (19°N) and American Samoa (14°S) [Conway et al., 1994, updated] delayed by 2 
months. The solid line is a linear fit to the long-term trend with a slope of 1.47 ppmv yr' 1 ; see 
text for details, (b) The derivative of the long-term trend (solid curve) shown along with the 
average growth rate calculated over the preceding 5 year (dot-dashed), 10 year (dotted), and 
15 year (dashed) periods. The horizontal line denotes the average growth rate from 1980 to 
1998 of 1.47 ppmv yr' 1 . 
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Figure 2. (a) The deseasonalized CO 2 boundary condition since 1880, derived for three peri- 
ods as described below. Period 1974-1999, from surface data for Mauna Loa, HI and Ameri- 
can Samoa, as described in the caption for Figure 1; period 1958-1974, based on annual mean 
CO 2 mixing ratios measured at Mauna Loa from 1958 to 1974 scaled to match the later years 
using overlapping data acquired from 1974 to 1993 [Keeling and Whorf, 1994]; period 1880- 
1958, using data from Antarctic Law Dome ice cores [Etheridge et al., 1996]. Symbols denote 
cores DE08 (triangles), DE08-2 (pluses), and DSS (solid circles). The dashed line is a linear fit 
to the boundary condition for 1980 to 1998 (same as the solid line in Figure la), (b) Test age 
spectra with mean age = 5 years. Spectra have the functional form of the analytic solution to 
the Green’s function for the 1-D advection-diffusion equation [Hall and Plumb , 1994; Andrews 
et al., 1999, 2001], which has been shown to fit age spectra from more complicated models of 
the stratosphere reasonably well, (c) CO 2 time series generated by convolving the boundary 
condition shown in Figure 2a plus a seasonally varying component (as in Figure 1 a) with the 
age spectra shown in Figure 2b. (d) CO 2 time series from Figure 2c detrended by subtracting 
the linear fit to the C0 2 boundary condition from 1980 to 1998 (dashed line in Figure 2a) de- 
layed by 5 years. This represents the differences between C0 2 evolution predicted using the 
long-term trend in Figure 2a and that which would occur if the CO 2 growth rate were constant 
at 1.47 ppmv yr' 1 with no seasonal cycle. 
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Figure 3. Average /lag (from (2)) versus N 2 O from the entire ER-2 data set, spanning lati- 
tudes from 70°S to the North Pole with data in all seasons except Southern Hemisphere sum- 
mer. N 2 O measurements are from the ATLAS instrument and have been normalized to the 
1997 tropospheric value following Strahan et al. [1999]. The solid circles represent the median 
value of /lag in overlapping 10 ppbv bins of N 2 O. The error bars represent 1 standard devia- 
tion. The solid curve is a polynomial fit excluding values corresponding to N 2 0 < 50 ppbv 
given by /Lag(N 2 0) = 0.0581(313 - N 2 0) - 0.000254(313 - N 2 0) 2 + 4.41x1 0 7 (3 1 3 - N 2 0) 3 , 
where 313 ppbv is the average tropospheric N 2 0 mixing ratio for the period from 
1992 to 1998. The open triangles represent mean ages for the midlatitude lower stratosphere 
from the empirical age spectrum analysis of Andrews et al. [2001], referred to as /spec in this 
paper. 
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Figure 4. (a) /lag (from (2)) versus N 2 0 for measurements from the NASA ER-2 aircraft 
within 2.5° of NASA Ames Research Center in Moffett Field, California (37.4°N), for each of 
the four seasons. N 2 0 mixing ratios have been normalized to the 1997 tropospheric value fol- 
lowing Strahan et al. [1999]. Each symbol represents the median mean age in a 10 ppbv N 2 0 
interval. DJF, December-January-February; MAM, March-April-May; JJA, June-July-August; 
SON, September-October-November. (b) The corresponding standard deviations, (c) Average 
/Lag as a function of N 2 0 for measurements from the NASA ER-2 aircraft within 2.5° of air- 
fields where the airplane was based during SPADE, ASHOE/MAESA, STRAT, and POLARIS 
and for measurements obtained within 5° of the equator. Each symbol represents the median 
mean age in a 10 ppbv N 2 0 interval, (d) The corresponding standard deviations. See text for 
discussion of sampling at each location. The curve shown in Figures 4a and 4c is the cubic fit 
to mean ages from the entire ER-2 data set, as shown in Figure 3. 
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Figure 5. (a) Mean ages from the tropical empirical age spectrum analysis of Andrews et al. 
[1999] for the seven tropical profiles obtained by the ER-2 from 1994 to 1997 plotted versus 
N 2 0 measurements from ATLAS. N 2 0 mixing ratios have been normalized to the 1997 tropo- 
spheric value following Strahan et al. [1999]. Shaded symbols denote tropical data obtained in 
March 1994 and February 1996. The solid symbols represent data from all other flights, which 
were obtained from NH summer through early winter (July - December). The solid line repre- 
sents the relationship between N 2 0 and mean age from the midlatitude empirical age spectrum 
analysis of Andrews et al. [2001]. The tropical age spectra give mean ages for March 1994 and 
February 1996 that are younger by as much as ~0.5 year than those obtained by applying the 
/spec:N 2 0 relationship derived for midlatitudes (dashed line), (b) N 2 0 profiles from ATLAS 
corresponding to the mean ages shown in Figure 5a. 
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Figure 6. (a) Variation with latitude of /"from all ER-2 observations of COt at 20 ±0.5 km 
from 1992 to 1998. Shaded symbols represent 10-s average data from each flight. Solid sym- 
bols with error bars (±2a) represent the average /"in overlapping 5° latitude bins, (b) Com- 
parison of average latitudinal variation of /"(calculated using (3), solid circles) and 71ag (cal- 
culated using (2), open triangles) at 20 ± 0.5 km from 1992 to 1998. (c) Seasonal variation of r 
as a function of latitude. 










Figure 7. Climatology of /"obtained by mapping the polynomial /TN 2 O relationship shown in 
Figure 3a onto the N 2 O climatology derived by Strahan et al. [1999] from in situ measure- 
ments of N 2 0 from the NASA ER-2 aircraft. Note that the N 2 O data are normalized to 1997 
levels as described in the section 2. 
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Figure 8. /Lag versus N 2 O for each of the (a) extratropical and (b) tropical OMS balloon 
flights; (c) average relationships between /lag and N 2 O for the extratropical and tropical 
flights. N 2 O data are from ALIAS-II, except for the November 11, 1997, flight where data are 
from the Argus instrument. N 2 0 mixing ratios have been normalized to the 1997 tropospheric 
value following Strahan et al [1999]. The polynomial fit to the ER-2 data from Figure 3 is 
shown for comparison, (d) /Lag versus CFC-1 1 for the OMS flights. 
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Figure 9. Scatterplot of mean ages derived from SF6 and CO 2 observations from the NASA 
ER-2. The dotted line is the 1:1 line, and the solid line is the 1.2:1 line. Shaded symbols rep- 
resent the 10-s flight data. Solid squares with error bars (±2a) denote the average SF6 age in 
overlapping 0.5 year intervals of CO 2 age. 
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Figure 10. (a) Comparison of CO 2 measurements from the OMS balloon flights with data 
from balloon-borne cryosamplers for NH midlatitudes. Midlatitude CO 2 profiles are typically 
characterized by a region of near-constant mixing ratio starting from 20-25 km and extending 
to the top of the profile (see, e.g., May 1998 profile in Plate 1). Symbols correspond to the av- 
erage mixing ratio in this region. The shaded line represents our stratospheric boundary condi- 
tion for CO 2 . The long-dashed line is a linear fit to the boundary condition from 1980 to 1998 
and the solid line has the same slope but with a delay of 4.5 years. Short-dashed lines corre- 
spond to delays of 4 and 5 years. Note that since the data here have not been corrected for CH4 
oxidation, the actual mean age is ~0.5 years older. See text for details, (b) Same data de- 
trended by subtracting the linear fit to the boundary condition delayed by 4.5 years (solid line in 
Figure 10a). 
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Figure 11. Average CO 2 mixing ratio corresponding to 215 ppbv N 2 O as measured on the ER- 
2 for 1992 - 1998 (asterisks), on the midlatitude OMS balloon flights, and as reported by Har- 
nisch et al. [1998] at ~20 km (simultaneous measurements of N 2 O were 215 + 5 ppbv); sym- 
bols for the balloon data are the same as in Figure 10. The shaded line and the dashed line are 
the same as in Figure 10. The solid line has the same slope as the linear fit to the boundary 
condition but with a delay of 3.5 years. 
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Plate 1 . Mean age profiles for each of the OMS balloon flights. The September 1996 and May 
1998 flights were from Fort Sumner, New Mexico (34.5°N); the June 1997 flight was from 
Fairbanks, Alaska (65°N); and the February and November 1997 flights were from Juazeiro do 
Norte, Brazil (7°S). 
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Plate 2. (a) CO 2 profile obtained on June 30, 1997. Magenta crosses represent air with “vor- 
tex-like” relationship between CH 4 and N 2 0. Data from the ER-2 for the same date are also 
shown. Approximate pressure altitude is calculated assuming a constant scale height of 7 km. 
(b) Relationship between C0 2 and N 2 0 for the same data and for the ER-2 flight of April 26, 
1997. (c) Relationship between f and N 2 0. 





